Lifetimes of agents under external stress 



OO 

o 

Q 

0^ 



I 

o 



> 
o 

^ ■ 
00 

ON : 

o ; 
•I— I 

^ : 

Oh- 

> : 
. 

X- 



Claus O. Wilke and Thomas Martinetz 
Institut fiir Neuroinformatik 
Ruhr-Universitdt Bochum 
(Submitted: December 8, 1998; Printed: February 2, 2008) 

An exact formula for the distribution of lifetimes in coherent-noise models and related models 
is derived. For certain stress distributions, this formula can be analytically evaluated and yields 
simple closed expressions. For those types of stress for which a closed expression is not available, 
a numerical evaluation can be done in a straightforward way. All results obtained are in perfect 
agreement with numerical experiments. The implications for the coherent-noise models' application 
to macroevolution are discussed. 

PACS numbers; 05.40.-a, 87.23. Kg 

Agents under externally imposed stress have been re- 
cently studied in coherent-noise and related models [Q-Q . 
These models display scale free distributions in a number 
of quantities, such as event sizes and lifetimes, or in the 
decay pattern of aftershocks. Coherent-noise models are 
very different from other models displaying scale free be- 
havior, such as sand pile models [||, as they do not rely 
on local interactions or feedback. Hence, they are not self 
organized critical. Considered the abundance of power- 
law distributed quantities in nature models such as 
the ones of the coherent-noise type can help understand- 
ing to what extent self-organized criticality is the right 
paradigm for describing driven systems, and to what ex- 
tent other mechanisms can provoke similar power-law dis- 
tributions. 

Despite the simplicity of the original coherent-noise 
model — agents have thresholds Xi ; if global stress exceeds 
a threshold x^, agent i gets replaced; with prob. /, an 
agent gets a new threshold — , no exact analytical results 
have been obtained so far. The distributions of event 
sizes and aftershocks have been studied in detail in ||] 
(event sizes) and in [0 (aftershocks). Both distributions 
can be regarded as being well understood. Nevertheless, 
the theoretical results are only of approximative charac- 
ter in both cases. 

In the case of the distribution of lifetimes, there are 
even less theoretical results. Sneppen and Newman ||] 
have given an expression based on their time-averaged 
approximation. This expression is right for certain 
stress distributions, as we will show below. However, it 
breaks down for slowly decaying distributions such as the 
Lorentzian distribution. Moreover, it is not clear when 
exactly it can be applied. 

In a recent paper , a different approach of calculat- 
ing the distribution of lifetimes has been taken, and the 
author claimed that the lifetimes obey multiscaling, with 
a, decrease for small lifetimes, and a. L^^ decrease 
for large lifetimes. Here, we will demonstrate that this 
statement is wrong. We will calculate the distribution 
of lifetimes exactly, without any approximations, and we 



will show that our results are in perfect agreement with 
numerical simulations. 

Our calculations are based on the observation that it is 
not necessary to know the distribution of thresholds p{x) 
for calculating the distribution of lifetimes. All we have 
to know is the distribution according to which agents 
enter the system, which is called pthrcsh(a;) in the nota- 
tion of and the stress distribution Pstrcss(a;). Once an 
agent has entered the system, it has a well defined life ex- 
pectancy, which is closely related to the probability that 
the agent will be hit by stress or mutation. Note that 
in this picture, we are considering only a single agent. 
Therefore, if we talk about lifetimes, it does not mat- 
ter whether the stress acts coherently on a large number 
of agents, or whether it is drawn for all agents indepen- 
dently. In this respect, the results we obtain in this work 
are of a much more general nature than the results found 
previously for event sizes or aftershocks. 

An agent with threshold x will survive stress and mu- 
tation in one time step with a probability p{x) equal to 



p{x) = (1 - /)[! -pmovc(a;)] 

= (!-/) / Pstrcss{x') dx' . 



(1) 



What is the distribution of the survival probabilities p7 
We denote the corresponding density function by u{p). 
Clearly, we have 



u{p)dp ^pt\ivcsh{x) dx 

— dx for < a; < 1 , 



(2) 



In the second step, we have assumed that the threshold 
distribution is uniform. This can always be achieved af- 
ter a suitable transformation of variables Hence, we 
find 



dx 

= Tp ■ 

The derivative dx/dp can be calculated from Eq. (|l|) 

dx _ 1 

dp (1 - f)Psti-css[x(p)] ' 



(3) 



(4) 
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and x{p) can be obtained from Eq. (|l|) by inversion. The 
density function is thus defined for p < Pmax, where 



Pmax = = (!-/) / Pstrcss(a;) dx (5) 
Jo 

stems from the condition that the thresholds are dis- 
tributed uniformly between and 1. Above Pmax, the 
density function u{p) is equal to zero. 

All agents with the same survival probability p gener- 
ate a distribution of lifetimes which reads 



g{L)=p^-Hl-p) 



(6) 



Here, g{L) is the probability density function for the life- 
times L. Note that the model works in discrete time 
steps, therefore the lifetimes L are integers, and g{L) is 
only defined for integral arguments. We can calculate 
the distribution of lifetimes h{L) by averaging over the 
distributions generated by different survival probabilities 
p, weighted with their density function u{p): 



h{L) 



u{p)p'^ \l-p)dp. 



(7) 



Equation (Q) can be explicitly evaluated for exponen- 
tially distributed stress, Pstrcss = exp(— a:;/cr)/CT. We find 



u{p) 



for < p < Pn 



with 



(l-/)[l-exp(-l/a)] 



(8) 



(9) 



After inserting this into Eq. (0) and doing some basic 
calculations, we obtain 



—j— + f<7 

^ -10 



dp . 



(10) 



It is possible to calculate the remaining integral with the 
aid of the identity (see [O, 15.3.1) 



t''-'{i-ty-''-\i^tz)-''dt = 



r(c) 



F(a,6; c; z) , 



(11) 



where F(a, b; c; z) is the hypergeometric function. We 
find 



h{L) = a 



-Pmax 



^ fIlA-.L+I:^"'^'' 



1-/ 



1-/ 



(12) 



The leading term ap^^^/ L is responsible for a, L ^ decay 
with cut off at L ~ 1/f. This behavior has been reported 



previously, and it corresponds to the approximation de- 
rived in [^. The correcting term vanishes with /. It is of 
importance only for extremely long lifetimes of the order 
1/f, for which it modifies the detailed cut off behavior. 

In Fig. 1^ we display Eq. ( p^ ) together with results from 
direct numerical simulations, for different values of /. 
The theoretical result is in perfect agreement with the 
measured distributions. The dependency of the cut off 
on / is clearly visible in Fig. |l]. 

Another stress distribution for which we can derive a 
closed analytic form for h{L) is the uniform distribution, 
Pstross(a;) = 1 for < a; < 1. We find 



u(p) 



1-/ 



for < 39 < 1 - / 



and 



fl_ f)L-i 



(13) 



(14) 



As in the case of Eq. (|1C|) , we get a leading and a correct- 
ing term. The leading term decays as with cut-off 
at L Ri 1//, and the second term modifies the cut-off 
behavior. Interestingly, the distribution of lifetimes is 
scale-free, although the distribution of event sizes in a 
coherent-noise model with uniform stresses is not a power 
law 1^. A plot of Eq. ( p^ is given in Fig. ||, together with 
the corresponding measured distribution. 

For the most other stress distributions, the integral in 
Eq. ^ can only be done numerically. This is the case, 
for example, for the Gaussian distribution, Pstress(2;) = 
^2/ (ttct^) exp[— a;^/(2(T^)]. Under Gaussian stress, an 
agent with threshold x will survive a single time step 
with probability 



pix) = il- f)cvi 



\V2cr 

where erf(a;) is the error function 
2 

erf(a;) = ^ / exp(-i^)rfi 
Vi" Jo 

Inversion of Eq. (|l5| ) yields 



(15) 



(16) 



(17) 



Here, by erf~^(z) we denote the inverse error function, 
obtained by solving the equation z = erf (a;) for x. We 
can calculate the density function of the survival proba- 
bilities with the aid of Eqs. (Q) and ([T7|). The resulting 
expression reads 



u{p) 



2 1-/ 



exp 



erf 



1-/ 



1 2^ 



(18) 



The numerical integration of u{p)p^ ^ (1 — p) is somewhat 
tricky for choices of a such, that Pmax/(1 — /) is very close 



2 



to 1, since the inverse error function has a singularity 
at 1. However, for moderately small cr, the integration 
can be carried out without too much trouble. The result- 
ing density function h{L) is shown in Fig. || for a = 0.15 
and / — 10"''. We find that, for L <C 1//, the function 
h{L) is almost linear in the log-log plot. A fit to the lin- 
ear region of h{L) gives an exponent r — 1.177 ± 0.01, 
which means h{L) decays slightly steeper than the 
decay predicted by the approximation of Sneppen and 
Newman. However, if we evaluate h{L) for much larger 
L and much smaller /, we find that the exponent r de- 
creases slowly towards the value 1 (Fig. ^). 

Let us now turn to the Lorentzian distribution 
Pstross(a;) — {2a/n)/{x^ + a?). In this case, a calcula- 
tion along the lines of Eqs. (|l|)-(||) yields the following 
distribution of survival probabilities: 



11 a 
HP) = TT- 



2 1-/ 



cos 



TT p 



21- f 



(19) 



Here, Pmax = (2/7r)(l — /) arctan(l/a). The result of the 
numerical integration is shown in Fig. ^. As in the pre- 
vious cases, we observe a perfect agreement between the 
analytic expression for h{L) and the distribution mea- 
sured in computer experiments. In the case of Lorentzian 
stresses, the distribution of lifetimes is clearly not scale 
invariant. 

In 1^] it has been claimed that the distribution of the 
agents' lifetimes under external stress decays as L^^ for 
small L. Among the four stress distributions considered 
in this work, we found a decay only for the uniform 
stress distribution. Hence the statement made in is 
wrong in general. We could verify the decay reported 
in 1^ for exponential or Gaussian stresses. As it was also 
stated there, the Lorentzian stress distribution does not 
produce a scale free distribution of lifetimes. 

A surprising result of this work is the observation 
that the properties of the distribution of lifetimes and of 
the distribution of event sizes in a coherent-noise model 
are largely independent from each other. We do find 
power-law distributed lifetimes under uniform stress, un- 
der which the distribution of event sizes is not scale free, 
and we do not find power-law distributed lifetimes under 
Lorentzian stress, which generates a scale free distribu- 
tion of event sizes. Consequently, we cannot infer from a 
power-law distribution of event sizes to one of lifetimes, 
and vice versa. Both distributions have to be investigated 
independently for every type of stress. 

Let us conclude with some remarks on the implications 
of our results for the application of cohcrcnt-noise or re- 
lated models to large scale evolution. In the context of 
macroevolution, the agents are regarded as species, or 
higher taxonomical units, such as genera or families [ p^ . 
The distribution of genus lifetimes in the fossil record 
follows either a power-law decrease with exponent near 
2, or an exponential decrease |l^,|l^. A decay can 



be observed in coherent-noise models with uniform stress. 
However, in this case the distribution of extinction events 
does not follow the decay - with s denoting the num- 
ber of families gone extinct in one time step - found in 
the fossil record j^] . The distribution of lifetimes closest 
to an exponential decay is, among the stress distributions 
we studied here, generated by Lorentzian stresses. But 
also in this case, the distribution of extinction events is 
significantly different from the needed decay of ex- 
tinction events. On the other hand, it seems to be typical 
for distributions generating a decay, such as expo- 
nential, Gaussian, or Poissonian, that the distribution of 
lifetimes decays as L~^. It is arguable whether any type 
of stress can actually generate the right type of distribu- 
tion for lifetimes and extinction events simultaneously. 
Hence, the coherent-noise models in their current for- 
mulation probably miss some important ingredient as a 
model of macroevolution. An effect which is not covered, 
and which has been shown recently to be of importance 
for the statistical patterns in the fossil record, is a de- 
cline in the extinction rate [p^-p^. For example, Sibani 
et al. |]l^,0 have demonstrated that the decay in 
lifetimes might be closely related to the decline in the 
extinction rate. 
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FIG. 1. The distribution of hfetimes of agents subjected to 
exponentially distributed stress, with cr = 0.08. The gray lines 
represent the results obtained from computer experiments, 
the black lines represent the theoretical prediction Eq. (p^. 
Theory and computer experiment are in perfect agreement. 
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FIG. 2. Comparison of theory and computer experiment 
for the uniform and for the Gaussian stress distribution. As 
in the case of Fig. |l|, we observe perfect agreement. 
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FIG. 3. The distribution of lifetimes of agents under Gaus- 
sian stress, with a = 0.15 and different values of /. For very 
small mutation rates /, the distribution of lifetimes becomes 
1/L for L larger than about 10^ time steps. 
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FIG. 4. With Lorentzian stresses, the distribution of life- 



times is no longer scale invariant. 



Ie2 le3 
time t (time steps) 



4 



